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We review a theoretical perspective of the dynamics of glass forming liquids and the glass tran- 
sition. It is a perspective we have developed with our collaborators during this decade. It is based 
upon the structure of trajectory space. This structure emerges from spatial correlations of dynamics 
that appear in disordered systems as they approach non-ergodic or jammed states. It is characterized 
in terms of dynamical heterogeneity, facilitation and excitation lines. These features are associated 
with a newly discovered class of non-equilibrium phase transitions. Equilibrium properties have 
little if anything to do with it. The broken symmetries of these transitions are obscure or absent in 
spatial structures, but they are vivid in space-time (i.e., trajectory space). In our view, the glass 
transition is an example of this class of transitions. The basic ideas and principles we review were 
originally developed through the analysis of idealized and abstract models. Nevertheless, the central 
ideas are easily illustrated with reference to molecular dynamics of more realistic atomistic models, 
and we use that illustrative approach here. 
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I. INTRODUCTION 

This article reviews current understanding of highly correlated dynamics characterizing structural glass forming 
liquids - super-cooled liquids approaching the glassy state. It is a hotly contested topic [H HI |3], and many recent 
reviews have been written about it, for example [U |5l [H [71 O |9l [10]. The majority of this writing adopts one of two 
perspectives: a thermodynamic view, the origins of which are often attributed to Goldstein [11] and to Adam & Gibbs 
[12] , and a small-fluctuation dynamical view called "mode coupling theory" (MCT) pioneered by Gotze [131 [TH [15] . 
Syntheses of these two approaches have been proposed [161 [13 [H] ^ some plausible arguments link aspects of such 
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syntheses to dynamical issues we discuss in this article [T9^. By and large, however, what we describe has a different 
basis, originating with discoveries that structural glass forming liquids exhibit pronounced dynamic heterogeneity. 

These discoveries were made with both experiment [20l EH [22l [21 [H [25] and computer simulation [26l [271 [28l [29l 
[30l[3l]. Dynamic heterogeneity emerges spontaneously from dynamics. It is a fluctuation dominated phenomenon 
that is largely independent of thermodynamics. Thermodynamics pertains mainly to mean behaviors. Mode coupling 
theory accounts for some fluctuation effects, but it is reliable only up to an onset of fluctuation dominance. As 
we describe in these pages, properties that distinguish structural glass formers from other materials are captured by 
physical models with large dynamic heterogeneity, irrespective of average thermal properties, and irrespective of many 
molecular details. 

Our presentation is qualitative, at the level of that found in the textbooks by Chandler j32j and by Barrat & Hansen 
[33] . We leave technical discussions of theoretical techniques and analyses to a future review [34]. Much of what we 
present here draws on results from numerical simulation of a simple atomistic model of a glass forming material. Our 
perspective on what occurs in that model and how it relates to experiment originated in theoretical work on idealized 
lattice models. These models are simple enough that they can be analyzed in detail, and remarkably, despite their 
simplicity, they successfully imitate many aspects of real super-cooled liquids, and many pertinent principles can be 
derived. 

Model building and solving is a traditional and important approach in statistical mechanics. Inarguable results for 
non-trivial models provides the basis for dismissing poor approximations, and for building the foundations of correct 
and general statements. For the topic of structural glass formers, the origins of the most used and useful models can 
be traced to Andersen and his co-workers [35l [36l [37] . One class is a collection of the so-called "kinetically constrained 
models" (KCMs) [35l [36l [38l [39l [40] , which are lattice models with Markovian dynamical rules that codify principles of 
facilitated [41] [42] and hierarchical [43] dynamics. The other class is composed of continuum models of fluid mixtures, 
which are atomistic models that can be simulated on computers [3T1|37]. At a superficial level, the two classes appear 
to be very different, but for their glassy behaviors they are similar. 

Approaching the end of this review, we reach the point where we are able to describe the recent discovery that 
models of structural glass formers exhibit a first-order phase transition between ergodic melt and non-ergodic glass 
phases [44 • It is a non-equilibrium phase transition, one controlled by fields that couple to trajectories - paths 
through state space. It is not controlled by thermodynamic fields like temperature and pressure, which couple to only 
states. Indeed, while only vague precursors are found in equilibrium measurements, a sharp unambiguous transition 
is found in systems driven away from equilibrium. This finding allows us to develop principles of the glass transition 
in terms of familiar concepts like order parameters, broken symmetries and wetting. 

II. DISORDERED CONDENSED MATTER 

Glasses are solids without evident structural order. They are formed from super-cooled liquids that slowly relax 
to equilibrium, sufficiently slowly that an arrested disordered material can be formed before an ordered crystal can 
assemble [4] |5] [6l [T] [HI [Sl [IH]- Figure [l] illustrates this juxtaposition of time scales. The transition to a glass is a 
dynamic process that occurs continuously over a range of temperatures, and the range depends upon the method by 
which the melt is cooled. The materials thus formed are truly solids, as much so as their ordered counterparts. 

Glasses do age and change their properties, some faster than others. Sometimes, this aging is the result of structural 
transformations, not fluid flow. For instance, initially clear kitchen glass wear may become foggy, and initially elastic 
organic glass may become brittle. These are examples of physical aging. The first is due to domains of quartz crystal 
growing within otherwise disordered regions. This phenomenon is accelerated through repeated heating cycles in a 
dishwasher. 

Solids, unlike liquids, are non-ergodic phases - atoms are fixed at or close to their initial positions for macroscopically 
long periods of time. In the case of crystals, the reason for this behavior seems obvious in view of the material's atomic 
structure. The energy and entropy are such that an ordered dense array of atoms is favorable, and an atom cannot 
move to a neighboring position without a massive rearrangement of atoms. The free energy for such a rearrangement 
would be formidable. Any appreciable atomic motions, therefore, must be due to defects that can move with lower 
free energy cost. Since a glass appears to be filled with structural defects, why is it also a solid? 

The answer that most often applies is appreciated by noting how the compression of a disordered system of hard 
space-filling objects may eventually reach a point of random close packing [47l [48]. This compression is sometimes 
referred to as "jamming" [49 . The resulting system is non-ergodic and therefore a glass. In a jammed or nearly 
jammed system, rare configurations that permit particles to move can represent only a small subset of deviations from 
crystal order. These rare configurations are the local dynamical excitations of a glass or glass former, excitations that 
lead to aging or relaxation. 

Behaviors of systems approaching the point of jamming can be observed in granular materials, and these behaviors 
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have commonalities with those of structural glass formers [501 ISH El]- In particular, restrictive forces between 
molecules in a super-cooled liquid lead to glassy behavior, and at high enough packing fractions the most important 
of these forces are repulsive intermolecular interactions. This idea, that the molecular structure and dynamics of dense 
disordered systems is often dominated by repulsive forces forms the basis for understanding the molecular behaviors 
of most normal liquids [53 , 54 , its origins can be traced to van der Waals [55^, to Bernal ^47j and to Widom [56^, and 
it underlies the standard equilibrium theory of simple liquids [571 EH] • 




FIG. 1: Temperatures and time scales for making glass. Panel (a) illustrates the time to nucleate and grow a crystal from its 
melt, Tnuc, and the time to relax the structure of a melt, r^. The former is a non-monotonic function of temperature T [45]. It 
is infinite at the melting temperature, Tm, and below that point it initially decreases before entering the regime where crystal 
growth is dominated by ripening; at that stage, its kinetics is limited by diffusion, and Tnuc then increases with decreasing 
T i46j. In contrast, Ta increases monotonically with decreasing T. The two times may intersect. The material will become a 
glass if the time scale over which the material can be observed and manipulated experimentally, texp, is shorter than Tnuc The 
glass transition temperature, Tg, is a temperature at which the structural relaxation time is about the same as the limiting 
experimental time. Tg is often taken to be the temperature at which the relaxation time Tq, becomes lO^s 31 151 IS]- is 
usually below the equilibrium melting temperature, Tm, and Tg depends upon both the system and the experimental time 
limits and protocols; for most observed glass formers Tg is about 2/3 of Tm [4l[5l[6]. Panel (b) illustrates typical temperature 
dependence of the enthalpy, if, for a liquid, its crystal and its glass. Curves of this type are constructed experimentally by 
recording the system's enthalpy as it is cooled or warmed at a rate of about l/texp- If a glass is formed, hysteresis occurs 
near T = Tg. The temperature at which the warming curve will have an inflection point is sometimes defined as the glass 
transition temperature, (c) Collapse of equilibrium transport data for fragile liquids in terms of the relaxation time at the 
onset to fluctuation dominance. To, the temperature at that onset. To, and a single energy parameter, J. The quantity r 
stands for either a structural relaxation time, t^, most often determined from dielectric measurements, or the shear viscosity, 
77, determined from transport measurements. The quantity To refers to this same property at the onset temperature. This is 
an application of Eq. ^ to more than 1000 data points from 67 different systems. Adapted from [68 . Inset shows the same 
data graphed in a so-called "Angell" plot [4j. 



III. REPULSIVE AND ATTRACTIVE FORCES IN CONDENSED MATTER 

Attractive branches of typical intermolecular potentials are responsible for making a dense liquid stable at standard 
pressures, but their effects are those of a mean or background potential that contribute little to energy fluctuations and 
forces that govern microscopic arrangements and movements of atoms [54l [56] . These movements and arrangements 
are much like those of a vibrating box of marbles or ball bearings. If the packing fraction is not too high, each atom 
can wiggle and jiggle and diffuse. Such random thermal motions are not rare (i.e., not activated), and they do not 



4 



involve correlated movements of many particles. At a high enough packing fraction, however, the space available to 
each particle is constricted, and any significant displacement of one atom must be accompanied by rearrangements 
of several neighboring atoms. Examples of these correlated motions are illustrated below when we discuss dynamic 
heterogeneity. 

Exceptions to where forces other than repulsions dominate are called "associated" liquids and "networked" glasses. 
Water is an associated liquid. Strong hydrogen bonds significantly infiuence its structure and dynamics. Most liquids 
do not possess such large attractive intermolecular forces. Hydrogen bonding causes liquid water molecules to arrange 
in a disordered tetrahedral network [59 . On small length scales, the liquid appears much like its ordered crystal ice. 
It is a relatively open structure. Entropic forces oppose it because it makes inefficient use of configuration space. 
On raising temperature, entropy quickly wins the competition with attractive hydrogen bonding. The preference for 
tetrahedral arrangements of water molecules ceases above about 50° C, well below the boiling temperature. 

Window glass is a networked glass. It is composed mostly of Si02 (i.e., silica), in which atoms interact with highly 
directional potentials that, like water, favor tetrahedral bonding of neighbors. Whether due to strong directional 
attractive forces, like those in silica, or due repulsive intermolecular forces that dominate in organic liquids, glassy 
dynamics and non-ergodicity are the results of molecular confinement - a constriction of accessible configuration 
space, and the concomitant sparsity of regions where molecular reorganizations are possible. This constriction is 
termed "frustration" [60 . In particular, part of the potential energy surface kinetically blocks or frustrates access 
to ordered lowest energy states, limiting configurations to higher energy disorganized configurations, of which there 
are many. In the case of organic glass, it is the irregular shapes of the mutually excluding molecules that create this 
frustration; in the case of networked glass, it is the directionality of attractive forces that create this frustration. 

While typical attractive forces in dense non-associated liquids cannot compete with repulsive packing forces, they 
can still have large effects. In particular, when molecular packing limits available configuration space to the point 
where structural reorganization involves coordinated displacements of several particles, many small attractive forces 
between separate pairs of molecules will act in the same direction and thus add significantly to the net forces controlling 
dynamics. In this way, attractive forces can enhance, not oppose, the effects of correlated particle motions. This is 
well known and understood in the context of equilibrium phenomena, where the presence of a large enough length 
scale heterogeneity, attractions are large contributors to the resulting unbalanced potentials between particles [61^ j62j . 
In glass forming liquids, the source of heterogeneity is correlated dynamics, which we will soon discuss in some detail. 
Nevertheless, in glass formers it is the repulsions that are the key frustrating forces, and the source of heterogeneity 
is correlated dynamics, which we will soon discuss in some detail. 

The terminology "packing fraction" refers to the fraction of space occupied by the molecular volumes of the particles 
in the system. For example, the packing fraction of an ordered array of close-packed hard spheres is (7r/6)pcr^ = 
7r/3\/2 ^ 0.7, where p is the number of spheres per unit volume, and a is the sphere diameter. Close packing of 
random hard spheres is estimated to be at about 10% lower density than that of ordered hard spheres An 
equilibrated ffuid of hard spheres freezes into its crystal when compressed to packing fraction of about 0.5, and the 
crystal's density is about 10% larger than that of the coexisting ffuid [64l[65]. Ordinary freezing of simple liquids is 
a reffection of this athermal ffrst-order phase transition of hard spheres [56l [66] . 

In real systems, repulsive forces are softer than hard cores, and as a result, the parameters characterizing molecular 
space fflling volumes are functions of temperature [54l El] • Roughly speaking, two atoms can approach each other only 
to the point where their repulsive potential of interaction reaches about /cbT (/cb stands for Boltzmann's constant). 
Thus, the typical distance of closest approach grows with decreasing temperature, so that at ffxed molecular density, 
p, the packing fraction increases with decreasing T. At constant pressure with decreasing temperature, there can be 
further increase in packing fraction because the density p will increase when the coefficient of thermal expansion is 
positive. 

Thus, it is clear from the context of intermolecular forces why super-cooling a liquid can produce constricting 
effects that slow dynamics and eventually trap it in a glassy state. But what is the nature of that transition? Figure [l] 
depicts what occurs with respect to changes in thermodynamic variables like temperature and pressure. At low enough 
temperatures (or high enough packing fractions), one observes hysteresis in thermal properties like the enthalpy. This 
behavior occurs smoothly, and it depends upon the rates at which the system is cooled and warmed (or compressed and 
expanded). It reffects intrinsic non-linear dynamics. For condensed matter, such non-linearity is generally associated 
with correlations between many degrees of freedom, and strong enough correlations produce order-disorder transitions 
like the sharp equilibrium transition between liquid and crystal. But the formation of glass exhibits nothing like this 
singular behavior, at least with respect to changes in thermodynamic variables like temperature and pressure. To 
reinforce this point, it is helpful to consider the experimental temperature dependence of structural relaxation, which 
we turn to now. 
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IV. GROWING TIME SCALES - UNIVERSAL, NOT SINGULAR AND NOT THERMODYNAMIC 

Structural glass forming liquids are generally catalogued as being either "strong" or "fragile." This terminology 
is due to Angell [4 . Strong materials are those for which measured equilibrium transport times grow with lowering 
temperature T in an Arrhenius fashion, 

log (r/To),t,„„g = (i?7fcB)(l/T - 1/To). (1) 

Here, we use r to stand for a structural relaxation time, Tq,, or a surrogate like viscosity, and Tq denotes a temperature 
at some reference point. When transport properties are well fitted by this form, provided E"" /k^T ^ 1, one may 
conclude that there is single type of event that leads to the longest time relaxation, and that the energy required to 
activate that event is ^* [32]. 

Fragile materials are those that are super- Arrhenius - their measurable relaxation times grow faster than exponential 
in 1/T. Most glass formers are fragile over the range of temperatures where they can be equilibrated (and thus observed 
to relax). The experimental data for these materials collapse to a form only slightly more complicated than Arrhenius, 
namely quadratic in reciprocal temperature [68^ , 

log (r/ro)fr,gi,e = {JiTof (To/T - if , (2) 

where is an energy determining the rate of growth of r as temperature is decreased from a reference temperature. 
To, which is called the "onset temperature." This temperature marks the onset to rapid growth of r upon lowering T. 
Transport properties of normal liquids have relatively insignificant temperature dependence, and this weak dependence 
mostly reflects the coefficient of thermal expansion, not activated dynamics [53", 54 . As such, the onset temperature 
is easily identified from data. Figure [TJc) illustrates the collapse of data below the onset temperature [68 . 

Temperature dependences of transport properties are frequently fit to the Vogel-Fulcher-Tamman (VFT) formula, 
log (r/rojyprp = A/{T — Tvft) 01 13 [6]. Here A, tq and Tvft are the three constants used to fit data to this 
form. These fitting constants are not directly related to properties of measurable reference points. For example, the 
temperature Tvft is the temperature at which r is imagined to diverge, hence the system cannot be observed to 
relax at that temperature. Despite this metaphysical quality, and despite the fact that experimental data make no 
compelling case for it [68l [76] , the VFT expression is sometimes referred to as a "law" . 

A. Dynamical perspective 

Super- Arrhenius behavior results from a distribution of relaxation times contributing to the longest time relaxation. 
For example, if the relaxation time roexp[A(l/T — I/Tq)] is random because the activation energy A is random, 
averaging over A can yield the quadratic function used to collapse data in Fig[l] In that case, A would obey Gaussian 
statistics with a mean of zero and a variance equal to 2J^ [69 . The super- Arrhenius behavior of Eq. Q also emerges 
from lattice models where microscopic reversible dynamics is hierarchical, with time scales for relaxing domains 
growing logarithmically with domain size [70 . Later, we will see hierarchical behavior in the dynamic heterogeneity 
of an atomistic model. 

As temperature is lowered, typical separations between mobile relaxing regions become larger, so that with hierar- 
chical dynamics, the number of coordinated steps required to relax an inactive region will also be larger. As a result, 
the apparent activation energy will grow with decreasing T. This growth will ultimately produce an activation energy 
that is larger than that required to simply avoid dynamical constraints, and in the absence of constraints dynamics 
cannot be hierarchical. Thus, when mobile regions are sufficiently sparse, i.e., when temperature is sufficiently low, 
we anticipate a crossover from super- Arrhenius to Arrhenius behavior, where the single pertinent activation barrier 
in the latter case is the energy required to disrupt the constraining network of forces. Lattice models predict this 
fragile-to-strong behavior [71, |72l [73]. Experimental observations of such a crossover have been reported for confined 
super-cooled water [74 , and for polymeric melts [75 , but the investigators attribute the crossover in the former 
case to a hypothesized equilibrium liquid-liquid phase transition, and in the latter case to supposed non-equilibrium 
phenomena. Clarification will require further experimental work. 

B. Thermodynamic perspective 

The VFT formula is often discussed in the context of the Adam-Gibbs expression, log {T/ro)j^Q oc 1/TScon-, where 
5'con is called "configurational entropy" [12]. 5'con is not uniquely defined, but it is usually equated with the molar 
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entropy of the super-cooled liquid less that of the crystal at the same temperature. This entropy difference can 
vanish, and the temperature at which it vanishes is Kauzmann's temperature, Tk [73 [78]. This would suggest that 
^con vanishes at T = Tk and that r would therefore diverge as T ^ Tk- Kauzmann's temperature would then be 
the temperature of an ideal glass transition - a non-zero temperature that would be unapproachable at equilibrium 
conditions. There is a convincing theoretical argument contradicting this metaphysical interpretation of Tk [79 . Yet, 
by adopting this approximation and expanding to linear order in T — Tk, we have 6'con ^ [{T — Tk)/Tk] AC, where 
AC is the difference between liquid and crystal molar heat capacities. The Adam-Gibbs expression then becomes 
the VFT formula, with A cx Tk/AC, and Tk = Tvft, thus suggesting a link between equilibrium thermodynamics, 
diverging time scales and the glass transition. Experimental support for this link is mixed [HOl [STJ [82l [83j , at least in 
part because applications require extrapolations of experimental data that exhibit hysteresis. 

Nevertheless, proposals relating the VFT formula and other divergent expressions to extrapolated thermodynamic 
properties are deeply entrenched in the literature on glassy physics. On theoretical grounds, this entrenchment may 
reflect results derived from mean field models of spin glasses [161 Ell ? which display both dynamical and thermody- 
namical singularities. But static disorder is a given rather than an emergent property of spin glasses, so a transition 
that might appear in spin glasses need not appear in structural glass formers. Not discounting the important practical 
role of contrasting and correlating thermodynamic and dynamic properties [85l |86l , there are too many unexplained 
exceptions to be able to use these correlations as the basis for a rigorous analysis. 

V. ORDER PARAMETERS AND TWO (OR MORE) STEP RELAXATION 

While the previous section emphasizes that there is nothing singular to be found in equilibrium data, there is in 
fact an order-disorder phenomena connected to the glass transition, but it occurs away from equilibrium, and it is 
not entirely controlled by thermodynamic variables like temperature and pressure. Explaining why requires some 
discussion, the first step of which is to consider appropriate order parameters. 

Just as the volume per molecule distinguishes two fluid phases - liquid and gas - there should be a characteristic 
property that distinguishes liquid and glass. A distinction based upon molecular configurations, however, is not 
obvious because liquids and glasses have similar if not identical densities and local structures. The most obvious 
difference between the two is dynamical. Normal and supercooled liquids are ergodic while glass is non-ergodic. An 
acceptable order parameter should highlight this distinction. For example, one measure of local dynamical activity is 
the squared displacement of a particle's position between time t and time t + At, the quantity |r^(t + At) — r^(t)p, where 
Ti{t) is the position of the ith particle at time t. At equilibrium, its mean value is independent of particle label and time 
origin, {\Ti{t + At) - Ti{t)\'^) = {\rj{At) - r^-(O)p), where iy^j. For large enough At, {\rj{At) - r^-(O)p) 2dDAt, 
where T) is a particle's self diffusion constant and d is dimensionality. The pointed brackets, (...), indicates equilibrium 
averaging over conditions at some point in time, or equivalently (assuming ergodicity) averaging over t for a very long 
trajectory. 

While a single-particle transport property like the self diffusion constant is a measure of dynamical activity, it is not 
by itself a suitable order parameter, which must be extensive in both space and time - extensive in space (or number 
of molecules) because we wish to distinguish systems viewed over macroscopic length scales (and thus detect broken 
spatial symmetry), and extensive in time because we wish to distinguish systems viewed over long observation times 
(and thus detect broken time symmetry). Such an order parameter can be formed by summing the local measure over 
the entire system for the entire time of observation (assuming both are very large). For example, one possible order 
parameter could be 

K = At^^\r,{t^At)-r,{t)\\ (3) 

t=0 1=1 

where N stands for the total number of particles considered, and tohs stands for the total amount of time the system 
is observed, incremented in the sum by steps of At. Another possible choice is 

Q = AtY,Y. ^^P{^^ • + - (^)] } ' (4) 

which is a measure of overlap of configuration space at time t with that at time t' = t + At, resolved on length scale 
27r//c. When 27r//c is a molecular length, extensive overlap for very large At is a signature of a non-ergodic solid [87] . 

For the case of equilibrium dynamics, the mean values of these two order parameters are given by standard equi- 
librium correlation functions. In particular. 



{K) = ATtobs {\ri{t + At) - ri(t)l') = 2dNtoi,sDAt, 



(5) 
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where the second equahty is true for At large enough for translational motion to be diffusive, and 

{Q)=Ntoi,,F{k,At), (6) 

where F{k^At) = (exp{ik • [ri(At) — ri(0)]}) is Van Move's so-cahed "self correlation function" or "intermediate 
scattering function" [[58j. 




FIG. 2: Equilibrium single-particle time-correlation functions for a two-dimensional glass forming mixture. The system is 
composed of 50% A-particles and 50% B-particles, which interact with WCA pair potentials [57|. The length parameter for 
the A-particles is the unit of length, a, while that for the B-particles is lAa. The pictured correlation functions refer to the 
larger of the two types of particles. (Those for the smaller particle behave similarly.) The net packing fraction of the mixture 
is Tvp (cr^ + a%) — 0.872. All particles have the same mass, m, and the same energy parameter, e. F(k, At) is the Van Hove 
correlation function defined in Eq. ([g]), and (Ar^(At)) is an abbreviation for a tagged particle's mean-square displacement 
(|ri(At) — ri(O)l^). Each of the upper two panels juxtapose behaviors above and below the onset to glassy dynamics. The 
bottom two panels juxtapose a measure of differing length scale dependences for dynamics above the onset temperature (bottom 
left) and below the onset temperature (bottom right). 



The terminology "equilibrium dynamics" refers to trajectories that are time reversal symmetric, preserve an equi- 
librium distribution of states, and, in the absence of external driving forces, reach this distribution from any initial 
condition. In this sense, a super-cooled liquid can be an equilibrium system. A glass, on the other hand, is the result 
of non-equilbrium dynamics, and we will see that its creation is the result of broken time- translation symmetry. 

Behaviors of the equilibrium correlations functions, and thus the equilibrium mean values of the order parameters, 
are illustrated in Fig. [2] The specific correlation functions shown are those computed by molecular dynamics simulation 
for a model d = 2 fluid mixture composed of classical particles interacting with the repulsive Weeks- Chandler- Andersen 
(WCA) potentials, i.e., the repulsive branches of Lennard- Jones potentials [54l[88]. Crystallization of this system is 
frustrated by the differing sizes of the two components in this mixture, one being 40% larger than the other. As such, 
it is possible to observe equilibrium dynamics in super-cooled fluid states of this system. The behaviors shown in Fig. 
[2] are typical of results obtained for any number of other such models, in dimensions two and three, that have been 
studied numerically [87^. 
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Figure [2] shows that a super-cooled glass former relaxes to equilibrium in stages. The very first stage is a fast 
inertial relaxation found in all fluids, even dilute gases. The next stage is the approach to a plateau that persists for 
more than an order of magnitude, at least as we see it at the lower temperature in the mean-square displacement 
and in At) for kcr = 27r. Notice, however, that this second stage is not evident at the lower temperature when 
coarse graining over length scales four time larger than <j, and it is not evident at the higher temperature. The 
wave- vector dependence of F{k^ At) therefore manifests some sort of heterogeneity in dynamics that is present at 
the lower temperature but not present at the higher temperature. The last stage of relaxation is the final decay of 
At). This final decay provides an estimate of a structural relaxation time, namely the 1/e-time for the Van Hove 
Function. We denote this time by re(A:), i.e., 1/e = F {k^re{k)); and we use simply Tg when k = 27r/cr, this value 
of k being roughly the location of the principal peak of the structure factor - the spatial Fourier transform of the 
equilibrium pair correlation function. 

While not graphed here in a way to highlight the functional from by which At) decays with time At, it is worth 
noting that at the lower temperature and k = 27r/(j, the final decay of the Van Hove function is a stretched exponential. 
But the stretching exponent depends upon length scale, so that for very small /c, At) decays with the normal 
exponential form characteristic of diffusive motion, exp(— /c^D At). This multifaceted behavior is a manifestation of 
fluctuation effects that we describe in the next section. 

The last stage of relaxation is called "a-relaxation." For a given glass former, structural relaxation occurs in a 
variety of ways, and there is a broad distribution of structural relaxation times. The wave- vector dependence of re{k) 
is one indication of this fact. The distribution of structural relaxation times can be usefully partitioned into at least 
two distinct classes of processes (see Fig. [s] and discussion surrounding it below). Irrespective of that analysis, there is 
a mean value for this distribution, and this mean is often termed the alpha relaxation time, Tq,. Tg can be a reasonable 
estimate of Tq,, but it is not the only possibility. Depending upon the physical issues at hand, Tg may or may not be 
the best estimate. 

The prior stage of relaxation, the second stage described above, is called "/^-relaxation" (or "fast" /^-relaxation to 
distinguish it from the slower "Johari-Goldstein" /^-processes observed in dielectric relaxation spectra [6]). The inertial 
and ^^-relaxations occur in both ergodic and non-ergodic phases — in crystals and glasses, as well as in equilibrated 
super-cooled structural glass formers. These relaxations reflect the residual dynamics of a rigid material. Unlike /3- 
relaxation, a-relaxation is absent from non-ergodic phases. In other words, non-ergodicity - broken time symmetry and 
rigidity for all times, what we describe as an "inactive" phase - is manifested by a finite value of (Q) /Ntohs = F{k^ At) 
for At ^ oo. This limiting value of the Van Hove correlation function is termed the "ergodicity parameter" or the 
"Edwards- Anderson order parameter" [87 . 

There are theories, most notably mode coupling theory (MOT) [T3[ HH [15] , that predict ergodicity breaking to occur 
at a non-trivial critical packing fraction or temperature as a consequence of non-linear feedback of large fluctuations 
in equilibrium dynamics. This dynamic singularity is not supported by experiment or simulation. Nevertheless, 
the location of critical packing fractions or densities and temperatures predicted from MOT serve as reasonable 
estimates for the boundaries between simple liquid behavior and fluctuation dominated dynamics [92j |93] , and the 
approach to that boundary from the simple-liquid side is well described by MOT [13^, |15]. It is a regime over 
which relaxation times change by several orders of magnitude, and MCT's successful description of this regime is a 
remarkable accomplishment of equilibrium molecular theory. 

VI. FLUCTUATIONS: DYNAMIC HETEROGENEITY, FACILITATION AND EXCITATION LINES 

Having discussed equilibrium mean values of order parameters, we now focus on fluctuations. This is the subject 
of dynamic heterogeneity. Spatial structures associated with fluctuations are resolved by relating parameters fields. 
For K and Q of Eqs. ([3| and Q, these are fields in space and time. For example, 

/> ^obs 

K= At V ^(r,t; At), (7) 

-^y t=o 

where integration extends over the net volume of the system, V, 

N 

/.(r, t; At) = ^ [r,(t + At) - r,(t)]' S (r - r,(t)) , (8) 

and S{t) is Dirac's delta-function, /^(r, t; At) is an example of what we call a "mobility" field. Its value at the space- 
time point (r,t) depends upon whether a particle is found at that point, and further whether that particle moves 
during the time period t to t + At. 
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Figure [3] illustrates this specific mobility field for different values of At in the dense equilibrium WCA mixture. 
(The initiation time t is unimportant because for large enough equilibrium systems, such as those pictured, any 
t is equivalent to any other t. In other words, equilibrium has time-translation symmetry.) Two temperatures are 
considered, one above the onset to glassy dynamics, and one below the onset. The most striking feature is the differing 
degrees of heterogeneity for the two cases. Panels at the same value of At/re illustrate fields at the same relative 
stages of average structural relaxation. The differing heterogeneities are thus not simply reflections of differing mean 
fields. 



At = O.lTe At = Te At = 5Te 




FIG. 3: Dynamic heterogeneity for the d = 2 WCA mixture considered in Fig. [2] The pictures are renderings of the mobility field 
/T:(r, 0; At) for typical equilibrium trajectories of 10^ particles. The displacement time, At is given in reference to the 1/e-time 
for the Van Hove self correlation functions at the corresponding temperatures. The rendering shades particles according to the 
size of the particle's displacement from its initial position. If the ith particle's displacement is nil, i.e.,|ri(At) — r^(0)| = 0, the 
particle is pictured as white (excepting a thin circle indicating its diameter). As the particle displacement grows, the particle 
acquires an increasing shade of grey, becoming completely black when \ri{At) — ri{0)\ ^ a. 



The darker the particles the more they have moved from their initial positions. After a long enough time At, all 
particles will move, and the pictures for both high temperature and low temperature cases will be uniformally black. 
But on the way to getting there, the low temperature mobility field is structured with thick connected mobile domains 
separating large immobile domains. The latter are cross sections of what we call "bubbles in space-time," about which 
we will soon say more. For now, notice that these large voids of immobility are present over extended periods of time 
in the low temperature case, whereas they are absent for all but the shortest At in the high temperature case. This 
juxtaposition indicates that spatial correlation of dynamics in the low temperature trajectory is much more significant 
than it is in the high temperature trajectory. 

"Large" is used here to describe lengths that are large compared to typical equilibrium correlation lengths, such as 
the coarse graining length over which the pair correlation function becomes structureless. This bulk correlation length 
is roughly 27r//Cmin, where /Cmin is the smallest wave- vector where the structure factor first deviates significantly from 
its limiting small wave- vector value. For dense fiuids (or glasses), this length is usually not larger than one or two 
diameters of the principal molecular species. In the high temperature system, one may infer from Fig. |3] that mobility 
coarse grained on this length scale spreads uniformly throughout the system so that its behavior is well approximated 
by the mean mobility. In contrast, in the low temperature system, motion does not spread so uniformly, and the 
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typical mobility coarse grained over the bulk correlation length is not well approximated by its mean. In other words, 
Fig. |3] demonstrates that dynamics of the high temperature system is like that of the mean field, while dynamics of 
the low temperature system is fluctuation dominated. 

Similar pictures of dynamic heterogeneity can be produced using any number of fields that capture local motion, 
for example, with the real part of Q(r, t; /c. At) with k ^ 27r/cr. This field is related to the net overlap Q in the same 
way that hz{T, t; At) is related to the net mobility K. Fluctuations in this field correspond to fiuctuations in the self- 
correlation function, and their structure factors, X4(<?5 ^7 ^t) = J dv' {Q(v^ t; /c. At) Q(r + r', t; /c, At))e*^ '^ , have been 
studied with computer simulations and theory by several workers jO^ 111 |95l [Ml EZl IMl IS3 fTOOl flOT] . X4(g; k, At) 
provides a useful quantitative measure of dynamic heterogeneity in an equilibrium system. It does not, however, 
provide a direct measure of space-time structure for transitions that break time-translational symmetry. For that 
case, the signature of a transition would be the divergence of mean square fiuctuations of extensive order parameters 
like K and Q, a divergence that grows super-extensively in tobs- From Eqs. ([3| and Q we see that this signature can 
result from correlations between fields like Q(r, t; /c. At) and Q(r', t^; /c. At) as \t — t'\ 00. The function X4^{Q'^ ^7 
contains information about t = t' only. 

A. Facilitated and hierarchical dynamics 

It is informative to enlarge Fig. [3] to the point where one can easily visualize the great extent to which particle 
are packed tightly. With that view in mind, it is interesting to know what is found from a movie of the mobility 
field with At growing continuously. Unedited, viewers are struck by high frequency motions, with particles fiipping 
rapidly back and forth between light and dark shades. These distracting high frequency motions can be filtered out 
by replacing Ti{t) and r^(t + At) in Eq. ^ with ri{t) and ri{t + At), respectively, where Ti{t) = {I /St) j^^ dtYiit). 
Rendering the resulting time-coarse grained mobility field, ^(r, t;At) with 5t ^ O.Olrg, proves illustrative. Movies 
created in this way [102 show the system possessed by low frequency low amplitude motions. These are the soft 
modes of a disordered solid. Darkening of particles signifying significant particle displacements takes place on time 
scales large compared to the apparent periods of those small amplitude motions, and this darkening takes place in 
the form of surges. String-like structures tens of molecules long emerge and then retract from central dark regions. 
The frequent retractions are indicative of the degree to which motion is restricted. An unconstrained random walk in 
two-dimensions would rarely lead to a perfect reversal over several atomic lengths. 

These surging strings are motions catalogued by Glotzer and her co-workers [30 . In a small system, these would be 
the only structures apparent before the entire system became dark. But for the size system pictured in Fig. [3) much 
more is seen. In particular, after several thrusts, a string emitted from one dark region no longer retracts, but rather 
blossoms into a larger fixed dark region that melds with the original dark region. From the new and now larger dark 
region, further strings surge, and the story repeats again at a larger length scale. These features exemplify a type of 
dynamics that is called "facilitated" and "hierarchical." It is facilitated because mobility in a region of space leads 
to motion (or relaxation) in an adjacent region of space. It is hierarchical because dynamics on smaller length scales 
(surging) is more frequent than motion on longer length scales (sticking and growing after several thrusts). Provided 
temperature is lower than the onset temperature, the great majority of events observed in renderings of ^(r, t; At) 
follow this type of path - facilitated and hierarchical dynamics. 

B. Excitation lines, decoupling and continuous-time random walks 

Another way to think about dynamical heterogeneity and facilitation is with the concept of excitation lines |103j . 
To the extent that dynamics is facilitated, excitations should form lines in space-time. This is true because facilitated 
dynamics requires an adjacent excitation for the birth of an excitation, and from time-reversal symmetry, an adjacent 
excitation is also required for the death of an excitation. Hence, space-time in a system with facilitated dynamics is 
structured with strings of excitations that are directed in time, and fiuctuations in numbers of excitations come from 
these lines dividing and coalescing. At any time slice, there is an equilibrium distribution of these excitations, and 
at equilibrium, this distribution is preserved from one time slice to the next. The dynamic heterogeneity of Fig. [3] 
should then emerge from overlaying the time series of these distributions. 

Figure [4] exhibits this type of structure. In this case, excitations are local regions of mobility, specifically the 
darkened regions located by the mobility field in upper-left panel of Fig. |4] These excitations are identified from 
the system's dynamics for the microscopic time 6t. For a glass forming liquid, they are sparse, like an ideal gas. 
When these dilute excitations at time t are added to those at times t + dt^ and those at time t + 2dt^ and so on 
up to those at t + nSt^ the resulting field is approximately the more expansive dynamic heterogeneity field shown in 
the lower left panel of Fig. [4j namely /^(r, t. At), where nSt = At. But one can also separately view each time slice 
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At = 10 5t 




FIG. 4: Excitation lines viewed from a bubble in space-time for the low temperature WCA mixture depicted in Fig. [3] Upper 
left panel shows the mobility field for the darkest (most mobile) particles when At = lOSt = O.lre. These dark regions are the 
elementary mobility excitations of the model. The right panel shows the trajectory of these excitations over a time span of 
t = 506t. The bottom left panel shows the mobility field when At = 506t which can be thought of as the accumulation of the 
mobilities from 100 times slices of the trajectory shown in the right frame. 



of the trajectory, /^(r, t,Jt) for t = Jt, 2(5t, n^t. This alternative produces the excitation lines of the super-cooled 
WCA mixture shown in the right panel of Fig. [4j On viewing it we can see that indeed the space-time structure this 
atomistic model of a glass forming liquid is as imagined from the perspective of facilitated dynamics. 

We first described this structure in Ref. |lQ3j , where we showed how its geometry provides explanations for quan- 
titative aspects of dynamic heterogeneity. In Ref. [103 , and in several subsequent papers, we have used kinetically 
constrained lattice models (KCMs), models where facilitation is a given or presumed property of dynamics. From 
Fig. |4j we see this property emerges from Newtonian dynamics of a sufficiently supercooled or compressed material. 
A few published papers have focused on how facilitation can be an emergent property jlQ4[ I1Q5( I1Q6[ I1Q7[ I1Q8( 1109] . 
The surging events described in the previous subsection coincide with the geometry of excitations in the KCM known 
as the "East" model [38 ; see, in particular. Figs. 2, 3 and 5 of Ref. [103^. It is therefore perhaps not surprising 
that the temperature dependence for relaxation times predicted by this simple model and its generalizations to higher 
dimensions [40 is the temperature dependence shown to agree so well with experiment in Fig. [l] 

One consequence of the excitation-line structure of space-time is that dynamical processes can be partitioned 
according to whether they coincide with exchange events or persistence events. The distinctions between these two 
classes of processes are responsible for a host of non-linear phenomena that are characteristic of deeply super-cooled 
liquids. To understand, consider a tagged molecule in a glass forming material. It might be a probe molecule that 
is investigated with single molecule spectroscopy, or it might be one of the molecules comprising the liquid itself. In 
view of Fig. (4) we see that over a time frame St it is most likely that this particle is immobile. In other words, it 
most likely sits in a void of mobility - what we call a bubble in space-time. Nevertheless, after some period of time, 
an excitation line will intersect an immobile particle, whereupon this particle will be able to move. It will continue 
to move until it steps out of the excitation line, or the line passes it by, whereupon the tagged particle will again 
be immobile. The picture we draw, therefore, is a dynamics where a tagged particle changes intermittently between 
mobile and immobile states [71 [23] . 
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FIG. 5: Probability distributions of log exchange time and log persistence time for the d = 3 WCA mixture for events defined 
as particles displacement of size a — 0.5cr. Adapted from Ref. ^91j where system details are described. The onset temperature 
for this system is To = 0.6e//cB [68]. As such, the left panel refers to a normal liquid temperature, and the right panel refers to 
a super-cooled liquid. 



The distributions of times characterizing this intermittent behavior can be viewed in terms of the time from its 
initial condition before it moves, and the time between two moves. The first of these is called a "persistence" time, tp, 
and the latter is called an "exchange" time, [HO . (In the literature on renewal processes, exchange and persistence 
times are sometimes called "waiting" and "excess life" times, respectively For uncorrelated sequences of events, 

persistence times and exchange times have the same distribution, and the probability that an event has not occurred 
in a time t is exp(— t/r), where r = (tx) = (tp). Such intermittencies are known as "Poisson" processes. In the case 
of glassy dynamics, however, event statistics is far from Poissonian because an event becomes more likely when a 
similar event has occurred more recently. In other words, facilitation implies that events are bunched in time. As a 
result, with the correlated dynamics of a super-cooled glass former, exchange times are typically much shorter than 
persistence times. Figure [5] demonstrates this behavior for the distributions of persistence and exchange times in a 
(i = 3 dimensional WCA mixture. At normal liquid conditions the two distributions coincide and are Poissonian. At 
supercooled conditions, however, the two distributions decouple with (tp) ^ {t^). 

One experimental manifestation of decoupling is the breakdown of Stokes-Einstein relations like Drj/T ^ constant 
or DTa ^ another constant. Here, D stands for a translational diffusion constant and 77 stands for shear viscosity. 
Because the average persistence time is the longer of the two mean relaxation times, it is an estimate of the structural 
relaxation time, Tq, ^ (tp). On the other hand, the time between random walk steps is so that D ^ o? /{t^)^ where 
a is a step size. As such, Dra ^ o?{t^) / {ty). The Stokes-Einstein relations hold when (tp)/(tx) ~ 1. This is the mean 
field result, and it is obeyed to a good approximation at standard liquid conditions. But in the fluctuation dominated 
regime of super-cooled liquids, the ratio can grow. The result [113' is a large failure of the Stokes-Einstein form that 
is seen experimentally |114[ I115| 1116] , and is illustrated for an atomistic simulation in Fig. [6] 

Typical excitation voids or bubbles of space-time persist for the time scale of structural relaxation, while particles 
diffuse by, in effect, surfing on the excitation lines that surround these voids. It is this surfing that makes D so much 
larger than anticipated by the Stokes-Einstein relationship. On the other hand, by applying a large enough external 
force to a tagged particle, the particle can be pulled out of an excitation line. In this way, a particle's drift velocity 
can actually decrease with increasing force. This negative response is yet another hallmark of fluctuation dominated 
correlated dynamics that has been observed in simulations. See for instance, Refs. |117( I118j . 

In cases of non-hierarchical dynamics (i.e., strong glass formers), fluctuation effects become irrelevant in dimensions 
d^2 [119J. For these systems, {t^) / {t^) will not grow upon lowering T. As a result, the predicted [73 crossover from 
fragile behavior to strong behavior upon lowering temperature might be tested by monitoring the behavior of Dtq^. 
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This possibility has been examined theoreticahy with KCMs |12Q[ 1121] , and the trends estabhshed with those models 
seem to be observed experimentally ^114[ I122j . 

There is a stochastic formalism that codifies our picture of the dynamics of a tagged molecule in terms of a series of 
periods of mobility or activity, with displacements larger than uninteresting small amplitude vibrations, punctuated 
by large quiescent periods of immobility or inactivity. It is the theory of continuous-time random walks (CTRW) 
|123L The formalism provides a means to distinguish the first from subsequent steps in a particle's trajectory. In 
the context of glassy dynamics, the first step coincides with persistence. It is the period before the tagged particle is 
intersected by an excitation line. The subsequent steps coincide with exchange events. 
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FIG. 6: Decoupling ratios for the d — 3 WCA mixture of Fig. [5] Adapted from Ref. [91]. We plot the quantity (tp)/(tx) 
(dark line) which measures deviations from Poissonian statistics a function of temperature. At high temperatures the two 
average timescales are the same and this quantity is unity. At low temperatures the persistence time rapidly decouples from 
the exchange time. We also plot the deviation from the Stokes-Einstein relation Dk^Ta, for k = 27r/(7. The decoupling between 
self-diffusion and alpha-relaxation timescale coincides with that of persistence and exchange. The onset temperature for this 
liquid mixture as it is defined by Eq[2] To = O.Ge/ZcB [BS]? is in harmony with the onset temperature for decoupling. 



The CTRW formalism has been applied to the Van Hove self correlation function with this perspective |124j , and 
it yields a result that is particularly vivid in the regime where (tp) (tx): 

F(k, t) ^ P{t) + [1 - P{t)] exp {-DkH). (9) 

Here, P{t) is the probability that a particle persists in its initial configuration for a time t, and D is the translational 
self-diffusion constant. The first term on the right-hand-side of ([9| is the persistence contribution to the mean 
concentration of overlap with initial conditions, the second term is the diffusive contribution. The balance between 
the two terms in ^ depends on the value of the wavevector k = |k|. That is, the relaxation timescale of self- 
correlations is lengthscale dependent |124i 

^e(fc)«(*p) + ^«(ip) + |^. (10) 



For large wavevectors, such as /c = 27r/cr in the WCA mixture, the first terms in rt9| and (10) dominate, so that 



Te{k) ^ Ta ^ {tp) . FoT Small cuough wavevector the second terms in d9b and (10) dominate, and relaxation becomes 
exponential with the diffusive (i.e., Fickian) relaxation time re{k) = 1/Dk^. Comparing terms of Eq. (10), we identify 



the characteristic lengthscale for the crossover between non-Fickian to Fickian regimes: /p = VDr^ oc \/ (^p)/(^x)- 
The lengthscale dependent timescales and the non-Fickian/Fickian crossover are central predictions arrived at from 
KCMs and facilitation theory, predictions that are in harmony with both experiments [125^ and simulations |126[ 1127] . 
The physical picture combined with the CTRW formalism has also has been used to analyze the X4-function [90], to 
interpret negative response [118^ , and to explain exponential tails [128^ observed in real-space self-correlation function 
^129 . It has also been extended to the study of "metabasin" transitions [130 . making a connection between KCMs 
and the energy landscape perspective [8 . 
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VII. ORDER-DISORDER IN SPACE-TIME 

The previous section described consequences of dynamic heterogeneity in structural glass forming liquids - broad 
distributions of length scales, time scales, decoupling, and negative response. Explicit pictures of dynamic hetero- 
geneity. Figs. [3| look much like those that one might see by studying the emergence or coarsening of an equilibrated 
state from a metastable phase [45]. It is as if there is an initial metastable immobile phase from which the stable 
ergodic or equilibrated phase emerges. In this section we discuss how indeed dynamic heterogeneity is the precursor 
to a phase-transition in trajectory space [1031 1138] . Doing so uses the so-called method of "large deviations" in 
generalizing traditional equilibrium statistical mechanics to a statistical mechanics of trajectory space 1131(1132] 1133] . 

We use the symbol Xt = (ri(t), r2(t), ...,ri(t), ...) to denote the point in phase space for the system at time t, and 
denote a trajectory with x{t) = (xq, x^t, ^25^, x^, x^^^J, where St is the time step for the trajectory, and the 
trajectory runs for a total of to\^^/5t time steps. For any system, there is an ensemble of possible trajectories. We 
use P[x(t)\ to stand for the distribution of trajectories at conditions of equilibrium. For deterministic trajectories, 
this distribution would be the equilibrium distribution of initial conditions times a product of delta-functions, one for 
each time step, specifying the rule by which the phase space point at time t is reached from that at time t — dt. For 
stochastic trajectories, P[x{t)] would contain a product of transition probabilities more general than delta-functions. 
(We are assuming dynamics is Markovian.) The specific form of P[x{t)] is not important. All that is important is 
that P[x{t)] is normalized, it preserves an equilibrium distribution, and it is time-reversal symmetric. 

The phenomenon we describe with this notation emerges from non-equilibrium dynamics - a non-ergodic glass 
phase that becomes the stable phase in an ensemble that is driven away from equilibrium. A way to construct 
non-equilibrium ensembles is to bias the equilibrium ensemble according to the value of some order parameter. For 
example, we can use the order parameter K[x{t)] given in Eq. Q and change P[x{t)] to 

Ps[x{t)] cc P[x{t)] exp {-sK[x{t)]) , (11) 

with the normalization constant 1/^s, where Zg = (exp {—sK[x{t)])). In analogy with equilibrium statistical mechan- 
ics, the parameter s is to K[x{t)] what 1//cb^ is to energy or what pressure times this factor is to volume. 

A positive value of s biases the ensemble towards trajectories that show little activity, a negative value of s biases 
the ensemble towards trajectories that show much activity. By adjusting the value of s one may push the ensemble 
far from regions that are most probable in the equilibrium distribution. In the limit of very large observation times, 
the quantity t~^^ In Zg converges to a time-independent function often termed the "large deviation" function. In other 
words, InZg is like a free energy of trajectory space, one that is extensive in time (and in space). A singularity in its 
density in the thermodynamic limit, [^"^^"bs ^s]v,tohs^oo^ such as a discontinuous derivative with respect to s, is 
the signature of a space-time phase transition with the chosen order parameter. Such behavior has been established 
to occur in several different KCMs, independent of any equilibrium phase transitions |134[ 1135] . 

The average order parameter in the non-equilibrium ensemble is 

Ks= -d\nZs/ds = Z-^ {K[x{t)]exp{-sK[x{t)])) . (12) 

In the context of molecular simulation, the second equality of Eq. ([T2| looks like a familiar expression used when 
re- weighting states in a thermodynamic perturbation theory calculation [32 l 1136] . The difference, of course, is that 
here we are averaging over trajectories, not states. As such, applications of these equations with numerical simulation 
require simulation techniques that do importance sampling of trajectory space. Transition path sampling [137 is 
a general methodology of this type, and it has been used in this way to implement the re- weighting of trajectory 
space |44 l 1138) Il39] . Figures [t] & [s] show results obtained for the Kob- Andersen mixture of Lennard- Jones particles. 

The evidence for a first-order phase transition induced by a small but positive value of 5 = 5* is compelling. In 
particular, on increasing system size and trajectory length, the change in Kg on passing from 5 < 5* to 5 > s* becomes 
more singular, and the order-parameter probability distribution becomes more bi-modal. Cleary, at 5 = 5* the system 
exhibits the signatures of macroscopic phase coexistence - a coexistence between an ergodic phase and a non-ergodic 
phase. At 5 = 5*, the material is not simply micro- heterogeneous, as it is at equilibrium. Rather here, pushed out of 
equilibrium, the domains of inactivity are as large as the total system. In terms of the equilibrium pair distribution 
function, ^(r), the structures of the equilibrium phase and the non-equilibrium phase are indistinguishable. Yet as 
measured by the van Hove correlation function, the non-equilbirium phase is non-ergodic - there is finite overlap with 
initial conditions for all times. 

For the KCMs that have been studied, the non-equilibrium active-inactive transition takes place at s* = 0+ |134( 
1135] , which means that equilibrium dynamics of these systems takes place at coexistence with the inactive dynamical 
phase. This location refiects that kinetic constraints in these KCMs are hard, i.e., they cannot be violated. For 
super-cooled liquids in general, kinetic constrains can be overcome, but in a way that becomes progressively unlikely 
as temperature decreases. One would therefore expect the coexistence point s* to be at a non-zero value, i.e. at 5 = 
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the inactive phase is destabihzed by constraint violating rearrangements. Nonetheless, at super-cooled conditions we 
expect 5* to be small, so that the equilibrium dynamics lies close to coexistence with the non-ergodic phase. 




FIG. 7: Evidence for first-order phase transition in space-time between an ergodic phase and a disordered non-ergodic phase. 
Average space-time order parameter Kg as a function of biasing field s, from molecular dynamics simulations of the Kob- 
Andersen Lennard- Jones liquid for N = 150 total particles and reduced density (Nact^ + NBcr%)/V = 1.05, at super-cooled 
conditions, T = 0.6e//cB. The onset temperature as defined in Eq. ^ for this system is To = O.Se/ZcB [68 . The field s couples to 
the net mobility K. As the length of trajectories increases, the crossover in Ks becomes sharper and happens at smaller values 
of s. For large tobs, the order parameter distribution at coexistence, s = s*, becomes bimodal, as expected for a first-order 
transition. While discriminating different regions of trajectory space with different values of s, the sampling avoids regions of 
trajectory space with ordered non-ergodic crystals by biasing against crystal order parameters. Adapted from Ref. [44^, where 
system size and other simulation details are described. 



Figures [7|& [8] show that the inactive states exist at temperatures below those of the onset of super-cooled behaviour. 
If one places the system in a configuration of the inactive phase, the system will remain inactive for times much longer 
than the typical structural relaxation time, even without biasing the dynamics. To the extent that s* goes to zero 
in the thermodynamic limit, the inactive phase can be called an "ideal glass". Its micro-states are vanishingly rare 
in the equilibrium ensemble, but once prepared, molecular dynamics will produce configurations that always overlap 
with that state. In other words, initial conditions are space-time surfaces which can bias the dynamics towards the 
active or inactive phase. But initial conditions taken at random from the equilibrium phase will essentially always 
select this equilibrium phase - the active dynamical phase [139 . 

To make this point explicit, imagine partitioning Zg according to its initial conditions, 

dxo p{xo) Zs{xo) = {Zs{xo)), (13) 

where p{xo) is the equilibrium distribution of initial states xq. As we have noted, in an ergodic phase, the net 
partition function, Zs is exponentially in trajectory length tobs- The same is true for Zs{xo). For large enough tobs, 



the annealed average over initial conditions, Eq. (13), is therefore dominated by those initial states with the largest 
values of (1/tobs) In ^^(xo). In this way, the non-ergodic phase can dominate even when its equilibrium weight is 
negligible. On the other hand, if the non-ergodic states have negligible equilibrium weight, the partition function for 
a quenched average over initial conditions, 

=exp[(lnZ,(xo))], (14) 

will always yield the partition function for the ergodic trajectories. In other words, the inactive phase can only be 
reached by some form of non-equilibrium driving |140j . 

This difference between quenched and annealed averages is a hallmark of rare region statistics [143\ Due to the 
proximity of bulk phase coexistence with the inactive phase, equilibrium trajectories display pronounced space-time 
bubbles of inactivity, and as demonstrated with KCMs, the largest or rarest of these inactive regions give rise to 
super- Arrhenius relaxation times and stretched time correlations [70l 1103] . Singular behaviors driven by surfaces 
imply surface phase-transitions [144J. For the transition between ergodic and non-ergodic phases, the "surface" is 
the initial conditions in space-time, and "wetting" is the overlap with initial conditions. From this perspective, the 
two-step correlators of Fig. |2] can be viewed as a pre- wetting profile, and the broken time-symmetry illustrated in Fig. 
[8] can be viewed as a wetted profile. 
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FIG. 8: Structural and dynamical measures of ergodic and non-ergodic phases of the d = 3 Lennard- Jones mixture of Kob and 
Andersen, for the same parameters as in Fig.[7| The field s couples to the net mobility K. (a) The radial distribution functions 
for the equilibrium ergodic phase, s = 0, and non-equilbrium non-ergodic phase, s > 0. (b) The Van Hove self correlation 
functions for the ergodic and non-ergodic phases, illustrated over a range of times very large compared to the equilibrium 
structural relaxation time Te. (c) The Van Hove self correlation function, illustrated over a range of the order of the equilibrium 
structural relaxation time, (d) Qualitative phase diagram depicting the ergodic-nonergodic coexistence line in the s — T plane. 
Adapted from Ref. [42, where system size and other simulation details are described. 



VIII. SUMMING UP 



The order-disorder phenomena discussed in the previous section is a bona fide phase transition between an ergodic 
state and non-ergodic state, with a partition function that is singular in a thermodynamic limit. It is, however, a phase 
transition that occurs away from equilibrium, with a partition function of trajectory space. With this transition in 
mind, it becomes comprehensible why a glass can form and why it may age very slowly. In particular, there are basins 
in trajectory space where initial conditions are remembered for all time. Those corresponding to ordered crystals are 
often accessible at equilibrium conditions. On the other hand, those associated with disordered solids are negligible 
at equilibrium but can become dominant in ensembles that are weighted away from equilibrium. Once driven to a 
configuration of one of the members of that ensemble, atomic configurations of that system will be remembered for 
times far into the future. These configurations are examples of non-ergodic glassy states. They can be metastable 
with respect to crystal phases, but reaching a crystal from one of these disordered states requires nucleation processes 
- creation of interfaces that grow to macroscopic sizes. The creation of interfaces are necessarily much slower than 
structural reorganization, and in the amorphous solid state, structural relaxation already takes a very long time. 

There should be no controversy surrounding the existence of this broken symmetry in time, a phase transition be- 
tween ergodic and non-ergodic phases. It is a fact established with exact analysis of idealized models, and convincing 
numerical evidence for more realistic atomistic models. The existence of this transition provides a natural explana- 
tion for the formation of glass and its equilibrium precursors, such as growing time scales, dynamic heterogeneity and 
intermittency, and it also explains why a glass and its melt are not easily distinguished by spatial structure alone. 
Nevertheless, it remains unresolved whether this transition is most often or even ever the actual glass transition 
encountered in real non-equilibrium processing of glasses. Resolution on this front must await the development of 
principles that link experimental protocols with those of re- weighting ensembles of trajectories, analogous to develop- 
ments of the last century that linked experimental protocols to Gibbs' statistical mechanics. We expect developments 
of this sort to begin to appear in the near future, possibly growing from the substantial and increasing body of 
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molecular theory and experiment on driven systems, such as viscous fluids under shear |145[ I146( I147j . 

A theme of this review is the recommendation to think dynamically not thermodynamically, or if the latter, to 
include time into the system's dimensionality. We flnd it remarkable that so many of the concepts of equilibrium phase 
transitions - order parameters, broken symmetries, wetting, and so forth - carry over to understanding order-disorder 
in space-time. The concepts we have described here for the highly fluctuating equilibrium dynamics of glass formers 
can also be applied to explicitly non-equilibrium problems such as that of aging systems jl48H149j . Furthermore, this 
understanding of dynamics may likely extend to areas beyond glassy physics, such as nanoscale assembly [15Q[ 1151] , 
protein folding dynamics and evolution [152 , and signal transduction in biology [153 . For example, signaling networks 
can be considered as many body problems of interacting proteins and protein complexes in the interior of the cell. 
While the collection of proteins and complexes is disordered and "crowded" [154j, the interacting units in individual 
signalling pathways are sparse and the interactions are often catalytic, i.e. facilitated, which is bound to give rise to 
fluctuation-dominated heterogeneous space-time structures analogous to those described above for the glass transition 
problem. This is an issue that has only has started to be addressed |155[ [T56j . 
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